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Abstract - We show how an appropriate stationary crystalline structure of the magnetic field 
can induce a partial fragmentation of the accretion disk, generating an axial jet composed of hot 
rising plasma twisted in a funnel-like structure by the rotation of the system. The most important 
feature of the obtained jet is its high degree of collimation, naturally arising from the condition for 
its existence. The presence of non-zero dissipative effects allows the plasma ejection throughout 
the axial jet and the predicted values of the accretion rate are in agreement with observations. 
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Introduction. — A relevant and puzzling feature of 
O compact astrophysical objects (e.g., Gamma Ray Burst [1] 
£jand Active Galactic Nuclei [2]) is their capability in gen- 
^erating highly energetic and collimated axial jets, well ob- 
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erved in the various electromagnetic bands. The way such 
jets remain collimated over very long astrophysical path is 
r ^ a non-trivial question, but indications exist that a pressure 
qq balance can take place as effect of the medium in which 
O^they are propagating. However, the mechanism which is 
l/^ responsible for the generation of such a peculiar emission 
of matter and radiation from a compact source remains al- 
t^most unidentified. A measurement of the difficulty to get 
<^> a satisfactory explanation of the jet phenomenon is pro- 
vided by the fact that one of the most promising proposal 
1 . I postulates the requirement of magnetic monopole effects 
>> [3-7]. Indeed, the possibility to formulate a reliable model 
k>for the generation of matter jet, without involving exotic 
j_j physics, must be regarded as a significant achievement in 
d understanding the behavior of the accreting material near 
a compact astrophysical source. 

A new framework for the investigation of the stellar 
wind and jet origin was offered by the studies pursued 
in [8,9], where it was shown how a local equilibrium con- 
figuration allows the existence of a periodic structure for 
the magnetic flux surfaces, as far as a proper account for 
the plasma magnetic backreaction is provided. This new 
approach, was extensively studied in [10-12] and general- 
ized to the global profile of the disk equilibrium in [13]. 
In this scheme, the equilibrium of a rotating stellar disk 



could favor the emergence of wind and jet seeds already 
in the steady state, as discussed in [14]. Such a work up- 
grades the original idea of an ideal purely rotating disk 
in which a crystalline magnetic field arises, by including 
non-zero poloidal velocities and matter fluxes. This fea- 
ture implies that the original local equilibrium must now 
deal with a non-zero azimuthal and electron-momentum- 
balance equations. The reason for which this scheme is 
able to account for high peaks of the vertical velocity (i.e., 
the seeds of winds and jets) consists just in the structure 
that the azimuthal component of the electron force bal- 
ance takes in the presence of poloidal velocities. In fact, 
since the azimuthal electric field must identically vanishes 
by virtue of the axial symmetry, one obtains a vertical ra- 
dial velocity that contains (in its denominator) the radial 
component of the magnetic field (almost vanishing in the 
background dipole field). Since the backreaction induces 
an oscillating radial profile in the radial component of this 
field (and, in the non-linear case, also in the disk mass 
density), it is immediate to qualitatively recognize that, 
in the O-point of the magnetic configuration, the vertical 
velocity takes extremely large values. This specific picture 
of the disk equilibrium configuration is certainly very in- 
triguing and promising, in view of the generation of winds 
or jets along the axial direction of an accreting compact 
object. Nonetheless, it contains three significant limita- 
tions: (i) the analysis is performed in a local radial sce- 
nario only, so preventing a localization of the peaks within 
the disk; (it) the emerging peak array appears to have a 
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periodic nature, in place of what is commonly observed in 
real sources; (Hi) the matter-flux lines are close, since they 
follow the magnetic profile and therefore no real ejection 
of matter is possible. 

The present analysis is just aimed at overcoming these 
difficulties and, by a global study of the plasma profile in 
the disk, we are able to construct an essentially single peak 
picture (the remaining ones being strongly suppressed out- 
ward). Furthermore, by including small dissipative effects 
in the plasma (according to the possibility to preserve the 
co-rotation theorem [12]), we show how the matter-flux 
lines become open allowing a real ejection of material out 
of the central object. In this respect, our work consti- 
tutes a significant step toward a definite model for a col- 
limated jet in the new scenario introduced by B. Coppi 
in [8]. Moreover, it offers a valuable theoretical frame- 
work for facing a more phenomenological characterization 
of the wind or jet generation in real astrophysical sources 
(in stellar, as well as, in galactic contexts). Regarding the 
phenomenological impact of the model, three main merits 
of this study must be put in order. First of all, the col- 
limated nature of the jet seed is guaranteed by the very 
short scale of the magnetic field crystalline structure in 
the disk (see [13]) and by the self-consistence of the equi- 
librium configuration equations. Since the region in which 
the radial component of the magnetic field (and also the 
mass density) is approximately zero is very tiny, then the 
outcoming jet is obliged to have a very small cross-depth. 
In the scenario of periodic jets, this fact is a shortcoming 
too, implying a small scale sequence of peaks and there- 
fore it is not clear how they appear when a macroscopical 
average is taken. Here the primary jet is only one and 
its small scale origin is directly related with its collima- 
tion. Secondly, it is possible to infer the location of the jet 
from the temperature profile of the disk. In fact, the lat- 
ter quantity is related to a parametric polynomial whose 
zeros fix the position of the jet. Finally, the present model 
is able to account for a non-zero accretion rate which is 
compatible with the typical observed values in real system. 

The paper is structured as follows. In the first Section, 
the basic equations and assumptions of the model for a 
thin accretion disk are introduced. In the second Section, 
the equilibrium equations are integrated and the behav- 
ior of the mass density and of the vertical velocity of the 
disk is outlined. As a result, vertical velocity peaks, cor- 
responding to the seed of winds and jets, are obtained. In 
the third Section, the consistency of the model is outlined 
and the temperature profile of the system is derived. In 
the fourth Section, we show how introducing dissipative 
effects (in particular, the resistivity) , the matter-flux lines 
becomes open implying an effective ejection of material. 
Brief concluding remarks follow. 

The Model. — In this Section, we construct a station- 
ary ideal magneto-hydrodynamical model for the equilib- 
rium of a thin accretion disk, surrounding a compact as- 
trophysical object of mass Ms- The disk is supposed to 



be embedded in the magnetic field B proper of the cen- 
tral object. In axial symmetry ((r, <f), z) being cylindrical 
coordinates), it writes 

B = (-d r ik/r; 0; d^/r) , (1) 

where ipQ denotes the magnetic flux surface associated to 
the background magnetic field. In what follows, the back- 
ground flux surface is taken in a reliable dipole form 

^ =Mor 2 (r 2 + z 2 )- 3 / 2 , (2) 

with p a = const. Moreover, we take into account that the 
accretion disk generates a backreaction magnetic field Bi, 
related to the corresponding magnetic flux function tpi: 

Bi = (-d r A/r; I/r; d z ^/r) , (3) 

where / = I{ip, z) defines the toroidal component and the 
total magnetic flux surface is described by tj) = ipo + 
with \ipi\ <C | Vo | - Let us now assume the following hi- 
erarchy in the plasma: the gradients of the flux-surface 
perturbations are of the same order or even greater than 
the background one and the second derivatives of V'l are 
always dominant. Considering the radial wave-number k 
associated to {i.e., d r ipi <~ kxpi), we moreover suppose 
to deal with a region of the accretion disk where kr 3> 1. 
In fact, it has been derived in [13] that, in the limit of thin 
disk (that we are considering), this region coincides with 
the whole disk. Our model is developed in the extreme 
non-linear regime for dipi 3> dipo- 

In this approximation scheme, the equations governing 
the dynamical system are the radial and vertical equilibria 
described by 

drP ~ 2pw ° ^ ^ + 1$ A ^ = ' (4a) 

+ =0, (4b) 

respectively (where A = d 2 + d 2 ), while the azimuthal 
equilibrium equation provides only the 0-component of 
the magnetic field and this feature is not needed for the 
addressed task. Here, p and p denotes the pressure and the 
density of the plasma, respectively, while lu represents the 
unperturbed angular velocity assumed to be Keplerian. In 
fact, since the gravitational field generated by the central 
object is described by the following Newtonian potential 
x(r, z 2 ) = GMs/Vr 2 + z 2 {G being the gravitational con- 
stant), the rotation writes 

, 2 = 4(r, Z 2 )EGM s (r 2 + z 2 r 3/2 . (5) 

As already discussed, inside the disk the magnetic field of 
the central object is well described the dipole- like config- 
uration Eq.(2). Thus, the following relation holds on the 
equatorial plane z = 0: 

u> 2 = c 2 (r, 0) - GMg/r 3 = GM S , (6) 
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where -00 = ipo{ r ity- Since the corotation theorem [15] 
states that the plasma angular velocity must be a function 
of the magnetic surfaces only, i.e., u = w(V>), it is natural 
to postulate 

u? = GM S V 3 /Mo . (?) 

extending Eq.(6) everywhere in the disk, also for a generic 
ijj (in view of the addressed perturbative approach in which 
l^i | -C \tpo |). It is worth noting that, comparing Eq.(7) 
and the form of u>k, the disk embedded in a dipole mag- 
netic field of the central object can not have a Keplerian 
behavior far from the equatorial plane. 

The pair of Eqs.(4) can not be solved analytically. In or- 
der to simplify the system, deriving Eq.(4a) with respect 
to z 2 and Eq.(4b) to r, in the limit of strong backreac- 
tion, one can obtain the following PDE by summing the 
resulting equations: 

2w — rd z 2 pVi) "o w d rP + ~. — 5 = , (8) 

where J*(-, •) is the Jacobian determinant, defined as 



J*(A,B) = det 



d r A 
d z 2A 



d r B 
d z 2B 



d r Ad z 2B - d z 2Ad r B , (9) 



(we note that, if A = A{B) than J*(A,B) = 0). In order 
to close the dynamical system, we now impose a constraint 
on tpi that, furthermore, allows to notably simplify Eq.(8). 
In particular, in view of the paradigm defined in [8, 9, 13] 
(according to which a crystalline structure of the magnetic 
field arises), we set 



A0i = A0i 



with A = const. < 



(10) 



that is commonly called the Helmoltz equation and we are 
allowed to neglect the last term in Eq.(8). A particular 
solution of Eq.(10) can be found as 



pi exp 



z 



sin(fcr) 
9{r) 
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where p± and h are constants, g(r) is a given function such 
that d r g(r) ~ g{r)/r (as we will see later, this function is 
connected to the temperature of the plasma), and where 
we have assumed that z 2 /h 2 <C 1. For this magnetic flux 
function, the eigenvalue is A = — (k 2 + 2/h 2 ). It is worth 
noting that a flux function of the form (11) generates, 
as expected, a magnetic field endowed with a crystalline 
structure. 

Considering ipi solution of Eq.(10), Eq.(8) reduces to 
a PDE for the function p only since J + = 0. In particu- 
lar, using Eq.(ll), the resulting equation for the density 
of the plasma can not still be solved analytically, due to 
the presence of the term d z 2 (pipi ) . In order to obtain an 
analytical solution, we can moreover assume 





Figure 1: Ideal plasma streamlines: #- function contour plot in 
the (x, it) plane where x = kr and u = z/h. 



where H denotes the vertical scale of the disk associated 
to the mass profile and h the vertical scale proper of the 
backreaction magnetic field. We can thus neglect the first 
or the second term on the right-hand side, depending if 
the backreaction is contained at all or is well outside the 
disk, respectively. 

The plasma momentum density pv obeys to the mass 
conservation equation V • (pv) = 0. Such a quantity can 
be thus expressed via a certain function as done for the 
magnetic field: 



pv 



-d z 0/r; pur; d r 0/r) 



(13) 



where 9 is called matter flux function. Writing the in- 
duction equation for the total magnetic field B(?/>) in the 
stationary limit, i.e.. <9 f B = V x (v x B) = 0, in terms of 
the functions 6 and ip, we obtain the following relation 



J(Mi) = 



(14) 



where we have implemented the extreme non-linear ap- 
proximation dipi ^> di/jg and the Jacobian J(-,-) is now 
defined as 



3(A, B) = d r Ad z B - d z Ad r B 



(15) 



s z2 (pV>iWi (p/H^ + p^/h 2 ) 



(12) 



From Eq.(14), we can get 9 — 9(ipi), i.e., the motion in the 
meridian plane takes place along the magnetic streamlines. 
The plasma is therefore driven by the magnetic field and it 
results to be frozen on the plasma streamlines preventing 
the plasma outflow from the system, as shown in Fig.l. 

The Emergence of a Jet. — Since the phenomenon 
we want to discuss is an axial jet of plasma extended out- 
side the accretion disk, and since the plasma is driven by 
the magnetic field only, we consider a vertical scale associ- 
ated to the backreaction magnetic field much greater than 
the one relative to the plasma density, i.e., h ^> H . 

In this approximation scheme, the first right-hand 
side term of Eq.(12) is leading and we can approximate 
d z 2(pipi) ~ ipid z 2p. We now take into account the thin 
nature of the disk H -C r, where H denotes, as already 
discussed, the vertical scale of the disk, i.e., the half-depth. 
Accordingly, the relation z<gr holds. In this respect, we 
can approximate wo — and Eq.(8) rewrites (considering 
Eq.(10)) 

d r p - 6r 2 0i d z 2p/u a = . (16) 
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Figure 2: Plot of the dimensionless density D of Eq.(17a) in 
the (x, u) plane where x = kr and u — z/h. 



Since the backreaction magnetic flux surface is separable 
in the form (11), the equation above admits the following 
solution 



p = p (r)D(r, z 2 ) = mr" d F(z 2 )N(r) , 



(17a) 



where po = mr d (with m = const, and d = const.) 
denotes the background mass density [13] and 



F = cxp 



2 /h 2 



6A£ 



N = exp 



kr 2 



cos(kr 



9{r) 



(17b) 



where (we remind that kr 3> 1) A = pi/po, and where 
£ is the so-called collimation parameter since it can be 
easly shown that v r /v z ~ £. This way, the more £ is 
small, the more the flux is collimated. Furthermore, from 
the expression of F, we obtain the half-thickness of the 
accretion disk to be defined as H ~ h^/6A£. Giving now 
a specific form to the function g, it results possible to 
determine the explicit form of the solutions. We infer a 
simple third-degree polynomial 



9(r) 



br 2 



(17c) 



with a, b, c constants. The choice of a polynomial is due 
to the request d r g ~ g/r and the third degree is the lower 
one that allows the presence of only one singularity in , 
placed in the origin of the coordinate system. 

The solution of the plasma density is therefore deter- 
mined and we now focus on the behavior of the vertical 
velocity: from Eqs.(13) and (14), it is easy to recognize 
that 

v z ~ (dr^/ipr) . (18) 

In Fig. 2 and Fig. 3, we plot the dimensionless plasma den- 
sity D = p/pQ (on the meridian plane) and the vertical ve- 
locity (on the equatorial plane), respectively. From Fig. 2, 
it can be shown how the crystalline structure of magnetic 
field induces a regular structure in density profile too. In 
particular, denser areas followed by emptier ones are ob- 
tained, i.e., a fragmentation. This behavior resembles that 
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Figure 3: Plot of the vertical velocity in the equatorial plane 
(z = u = 0) in function of x — kr (we have set 9 = (C/k)ipi 
with C = const.). 



seen in [8]. In Fig. 3, as a result, we find a space region 
(in correspondence of the density minima) in which the 
plasma develops a very large vertical velocity that repre- 
sents the formation of a jet. The jet is carried outside the 
plasma by the magnetic field reaching a height h but the 
motion results closed. In fact, due to the magnetic field 
action, the plasma is forced to move toward the equato- 
rial plane as shown in Fig.l. Of course, this feature is in 
contrast with observations and we will fix it in the last 
Section. Finally, in Fig. 4 is depicted a three-dimensional 
view of two jet streamlines, twisted together by the ro- 
tation of the accretion disk. It is shown the funnel-like 
nature of such an axial jet and the high collimation of the 
stream. 
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Figure 4: Three-dimensional picture of the jet, in which only 
two different plasma trajectories are represented. 



Thermodynamical Features of the Jet. Let us 

now focus on the thermodynamical properties of the sys- 
tem. In particular, substituting the expressions of p and 
rpi (given by Eqs.(17a) and (11)), in the system of the 
radial and vertical equilibrium Eqs.(4), one gets a couple 
of PDEs for the function p only. Since we deal with two 
differential equation for just one unknown function, the 
compatibility of such an overdetermined system has to be 
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discussed. First of all, we define the following dimension- 
less quantities 

x = kr , u = z/h , Y = kipi/d r ipo , (19a) 

K = GMsm/fio , A = {Khy 2 , P = p/B 2 (z) , (19b) 

P = 2A(P + (2£ - 1)Y 2 / 871 ") • (19c) 

Setting d = 3 in the definition of po, the system (4) 
rewrites now as 

d x P-6£DY = 0, d u 2P + D = 0. (20) 

From the second of Eqs.(20), we obtain the condition 

P = -N(x) jF{u 2 )du 2 , (21) 

and, substituting the expression above into the first of 
Eqs.(20), the compatibility condition results simply £ -C 1, 
implying a well collimated flux. We underline how this fea- 
ture is crucial, since it has been obtained just by requiring 
the model consistency. It is a self-consistent behavior of 
the system. 

Let us now discuss the morphology of the plasma pres- 
sure obtained from P, and, accordingly, of the temper- 
ature T of the accretion disk. In fact, these quantities 
are connected by the perfect gas equation of state, fixed 
as p = v 2 p ~ pT, where v s indicates the plasma sound 
velocity. In Fig. 5, we plot the plasma dimensionless tem- 
perature field T = P/D. It is easy to realize that, in 
correspondence to the velocity peaks (or the density min- 
ima) , temperature peaks are outlined as well. This implies 
that the jet seed is characterized by an upward motion of 
very hot plasma. In this respect, the temperature field can 




Figure 5: Plot of T = P/D on the (x, u) plane. The parameters 
are set as for the previous analyses. 

be regarded as responsible for the presence of the velocity 
peak itself, i.e., for the particular form of the g-function 
introduced in Eq.(ll). If the disk possesses a region in 
which the plasma is hotter with respect to the mean tem- 
perature, a strong modulation of the crystalline structure 
appears, generating the ring decomposition and the ver- 
tical jet. Moreover, the presence of a hotter part in the 
plasma is clearly more probable in the inner regions of the 



disk and this results to be in agreement with the obser- 
vations. The temperature profile grows with the vertical 
coordinate and this is due to the fact that the jet is rarefied 
in the outer part of the system and because no cooling-by- 
radiation mechanism has been introduced in the model. 

Dissipative Effects. — In the first Section, it has 
been outlined how the plasma velocity (related to ff) is 
function of tp and it results to be parallel to the magnetic 
field. The motion of the jet seed is therefore closed since 
the particle trajectories are frozen on the closed magnetic 
lines. In this respect, in order to obtain a jet implying an 
effective matter transport outward the disk, a mechanism 
responsible for the misalignment of the gradient of the flux 
functions 9 and tp is requested. To this end, we introduce 
in the equilibrium equation a viscous and a resistive coef- 
ficients rj v and r) r , respectively. We require these effects to 
be relevant in the magnetic field induction equation only. 
This is a regime in which resistive effects are stronger than 
viscous ones (also corresponding to the limit of a magnetic 
Prandtl number greater than one). Moreover, we under- 
line that in [12] it has been shown how the corotation 
theorem [15] holds under the hypothesis of stationarity up 
to the second order in small dissipative effects, poloidal 
velocity and toroidal component of the magnetic field. In 
our approximation scheme, one obtains the following mag- 
netic induction equation 

V x (v x B + ?7 r V x B) = , (22) 

while the system (4) remains unchanged. The vertical and 
radial components of the equation above can be combined 
to obtain, in the extreme non-linear regime we treated, 

J(ip 1 ,6) + pr ]r rAij 1 =0, (23) 

in place of Eq.(14). It is important to stress that the pres- 
ence of a non-zero constant resistivity coefficient implies 
the misalignment discussed above. In fact, we obtain the 
relation \Aipi x A6\ = —r) r rpAipi, that guarantees the 
emergence of open flux lines of matter. Eq.(23) corre- 
sponds to a PDE for the ^-function only, since the radial 
and vertical equilibrium equations remain unchanged (be- 
cause of the smallness of the poloidal velocities) and we 
can simply import, in the dissipative scenario, the func- 
tions ipi and p obtained in the ideal case. Using now 

d(r ) z) = 6 Q (il>)+9 1 (r,z), Q = —y/lh, (24) 

XT] r p 

(where 9± must be an odd z-function in order to obtain a 
non-zero accretion rate of the disk) and the dimensionless 
functions (19), Eq.(23) rewrites 

J(Y, 9) = (1 - 2£)Y , (25) 

where J denotes the Jacobian determinant expressed in 
dimensionless coordinates. Eq.(25) is an equation for 8 
only and can be solved numerically yielding the plasma 
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streamlines plotted in Fig. 6. As discussed above, the pres- 
ence of an even small constant resistivity coefficient opens 
the plasma streamlines allowing an effective plasma out- 
flow through the jet. The opening of the streamlines is 
shown in correspondence of the velocity peak of Fig. 3. 




Figure 6: Dissipative plasma streamlines: ^-function contour 
plot in the (x, it) plane where x = kr and u = z/h. 

The integrated 0-function allows us to calculate the ac- 
cretion rate towards the central object since the following 
relation holds 

M(r) = —2nr / pv r dz ~ — 2irr pv r dz — 4ir8i(r, h) . 
J-H J-h 

Considering a stellar astrophysical scenario, at the bound- 
ary layer between the disk and the star we thus obtain 

M ~ 1(T 1O M /Year , (26) 

a value compatible with observations [3]. In this estima- 
tion, we have considered a resistivity coefficient of the form 
T] r ~ I0 12 (T/K) -3 / 2 cm -2 s _1 [16] and values for mass, 
radius and temperature typical of the accretion disks as- 
sociated to a cataclysmic variable, i.e., A/disk = 1O _4 M0, 
i? disk = 10 6 Km and T d isk = 10 7 K. 

Summarizing, we have developed a model of a station- 
ary jet phenomenologically viable, at least on a qualita- 
tive level. In fact, by virtue of the dissipative effects, the 
velocity field parts from the magnetic field describing a 
self-collimated flux directed outwards the system and the 
accretion rate associated to this feature is compatible with 
observations. 

Conclusions. — As in [13] the crystalline structure 
of the magnetic field in a thin accretion disk (reacting to 
the dipole field of the central object) has been extended 
from a local to a global profile, also in this paper we have 
upgraded the analysis of wind and jet seeds performed 
in [14] toward a global picture. This allows to reproduce 
important observational features like the isolated nature of 
the jet and a non-zero accretion rate of the disk as a whole. 
Within such a theoretical paradigm, we have developed a 
stationary equilibrium configuration of the disk which is 
able to incorporate the presence of an highly collimated 
jet profile, corresponding to those region of the disk where 
the matter density acquires an absolute minimum in its 
crystalline radial dependence. 

The main assumptions we have dealt with are the 
small values of the poloidal velocity and of the requested 



Prandtl number, according to the preservation of the co- 
rotation theorem. A non-zero resistivity coefficient has 
been needed to account for real ejection of material from 
the disk, since the matter flux lines are no longer isomor- 
phic to the magnetic field profile. The model contains 
also a free radial polynomial function g(r). which fixes the 
position of the jet (g = 1 would imply that the jet be lo- 
cated too close to the central object, i.e., almost at r = 0). 
However, such a function can be directly connected to the 
thermodynamical properties of the region where the jet 
arises, in particular with its temperature. The radial do- 
main where the jet takes place corresponds to a peak of 
the temperature profile, which, in turns, increases with 
the distance from the equatorial plane. 

The most appealing feature of the proposed model is 
that the jet collimation comes out from the compatibil- 
ity of the equilibrium configuration system, and, at the 
same time, the jet radial extension is very tiny due to the 
small scale of the crystalline profile. This feature, associ- 
ated to a non-zero accretion rate of the disk (the values 
are in the observed range as soon as the typical orders of 
magnitude of a stellar system are considered), makes the 
jet model here derived of very promising phenomenolog- 
ical impact and surely an intriguing perspective in view 
of a non-stationary extension of the equilibrium config- 
uration. The main merit of this work is the significant 
upgrading provided with respect to the local analysis [14] , 
which enforces the idea that the crystalline profile of the 
backreaction is a suitable scenario to implement the jet 
formation from an axisymmetric accretion structure. 
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